




ERRATA NO. 1 


NACA REPORT 1028 

EFFECT OF ASPECT RATIO ON THE AIR FORCES AND 
MOMENTS OF HARMONICALLY OSCILLATING THIN 
RECTANGULAR WINGS IN SUPERSONIC 
POTENTIAL FLOW 
By Charles E. Watkins 

1951 


Page 7 , column 1: In equation ( 38 a) , the last bracketed expression 

should be corrected as follows: 

[3p 2 + 6 - Ax q (2 + p 2 )] 

Page 7 , column 2; In equation ( 38 b), the last bracketed expression 
should be corrected as follows: 

|7p 2 + ^ - 5xq)J 

Page 8 , column 2: In equation (43a) the last term should be corrected 

as follows: 

}f > s + 2 s 

Page 9 , column 2: In equation (45c) , the factor 2 preceding the second 

parenthesis should be deleted; that is, the second term within the 
bracket should read 

- (1 - 2xo)(f x + G X ) 


Page 9, column 2: Equation (46a) should be corrected to read as follows: 



It is pointed out that the foregoing errors have been corrected in a 

subsequent NACA publication (NACA TN 3076 by Nelson, Rainey, and Watkins), 
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LAMINAR BOUNDARY LAYER ON CONE IN SUPERSONIC FLOW AT LARGE ANGLE 

OF ATTACK 1 

By Franklin K. Moore 


SUMMARY 

The laminar boundary -layer flow about a circular cone at 
large angles of attack to a supersonic stream has been analyzed 
in the plane of symmetry by a method applicable in general to 
the flow about conical bodies . 

At the bottom of the cone , velocity profiles were obtained show- 
ing the expected tendency of the boundary layer to become thinner 
on the under side of the cone as the angle of attack is increased . 

At the top of the cone , the analysis failed to yield unique 
solutions , except for small angles of attack. Beyond a certain 
critical angle of attack , boundary -layer flow does not exist in 
the plane of symmetry , thus indicating separation . This 
critical angle is presented as a function of Mach number and 
cone vertex angle . 

INTRODUCTION 

The supersonic aerodynamics of pointed bodies has con- 
siderable current interest in connection with the design of 
aircraft and missile fuselages. An important feature of the 
flow about such bodies is the behavior of the boundary layer 
and, in particular, the flow separation which may occur along 
the low-pressure side of the body due to angle of attack. 
The present report will consider the development of the lami- 
nar boundary layer on the surface of a right circular cone at 
an angle of attack to a supersonic stream (see fig. 1). The 
conical configuration may be considered an idealization of 
the nose portion of a supersonic aircraft fuselage. 

Outside a thin boundary layer on a cone, the nonviscous 
supersonic flow (upon which the boundary layer itself de- 
pends) is “conical” in the sense that physical quantities 


(such as velocity and pressure) are constant along any ray 
proceeding from the cone apex. The description of this 
outer flow, contained in references 1 to 3 and elucidated in 
reference 4, is considered adequate for the purposes of this 
report, but subject to restrictions which will be discussed 
subsequently. 

In figure 2 is shown qualitatively the circumferential 
pressure distribution on the cone surface predicted for var- 
ious angles of attack (see ref. 4). These pressure distribu- 
tions depend only on the character of the nonviscous flow 
beyond the boundary layer, on the assumption that the 
boundary layer is extremely thin. When the angle of attack 
is very small, the pressure decreases monotonically from 
the bottom of the cone around to the top. For larger 
angles of attack there appears a region near the top of the 
cone wherein the pressure gradient reverses and the pressure 
increases toward the top. As the angle of attack is further 
increased, this region becomes greater in extent. 




(a) Small «. (b) Moderate <*. (e) Large a. 

Figure 2. — Pressure distribution around cone for various ranges of 
angle of attack. 



As a consequence of the conical nature of the nonviscous 
flow, it is shown in references 5 and 6 that the laminar 
boundary layer has parabolic similarity along generators of 
the cone; that is, velocity, pressure, and density inside the 
boundary layer are constant along any parabola (see fig. 1) 
drawn in any one meridional plane (plane passing through 
the body axis) . Of course, circumferential variation of these 
quantities is to be expected when the cone is at angle of 
attack. 

In reference 7 the effect of angle of attack on the laminar 
boundary layer is analyzed, in the limit of very small angle 
of attack, with the result that the boundary layer tends to be 
thicker on the top of the cone than on the bottom (fig. 3(a)). 


i Supersedes NACA TN 2844, “Laminar Boundary Layer on Cone in Supersonic Flow at Large Angle of Attack” by Franklin K. Moore, 1952. 


1 



2 


REPOET 1132 — -NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



(o) (b) (c) 

(a) Small a. (b) Moderate a. (c) Large a. 

Figure 3. — Cross-sectional views of boundary layer on cone at various 
angles of attack. 

This is to be expected since the fluid near the base of the 
boundary layer has low inertia, is therefore inclined to 
follow the direction of the circumferential pressure gradient 
more closely than is the outer flow, and thus tends to drain 
away from beneath the cone and accumulate near the top. 
No separation is encountered because, for small angle of 
attack, the pressure gradient is always favorable (fig. 2(a)). 

For larger angles of attack, when the pressure gradient 
reverses direction near the top of the cone (fig. 2(b)), experi- 
ment indicates the formation of boundary-layer “ lobes’ ’ 



(fig. 3(b)). When the angle of attack is further increased, 
the lobe pattern finally breaks away from the body to form 
a vortex street (fig. 3(c)). 

Recently (ref. 8), at the Lewis laboratory a brief experi- 
ment was carried out in which a total-head probe was placed 
near the surface at the top of a cone and pointed toward 
the cone apex (fig. 4(a)). The cone was mounted in a 
supersonic wind tunnel, and the probe was used to measure 
the total head in the boundary layer at a fixed height above 
the surface as the angle of attack was varied by rotating 
the cone in the meridional plane containing the probe. 
Figure 4(b) shows the result of this test. The decrease in 
indicated total head as the angle of attack was increased 
from negative to positive values may be interpreted to 
mean an increase in boundary-layer thickness at the top 
of the cone. Beyond a certain angle of attack, this tend- 
ency reverses, and the boundary layer apparently becomes 
thinner as the angle of attack is increased. This is a pos- 
sible indication of the tendency to form lobes, as illustrated 
in figure 3(b). 

In the present report of research conducted during the 
summer of 1952 at the NACA Lewis laboratory, the laminar 
boundary layer in the meridional plane of symmetry of the 
flow is analyzed for large angles of attack in order to provide 
velocity profiles on the bottom of the cone and to provide 
a certain degree of insight into the question of separation 
on the top. 

SOLUTION OF BOUNDARY-LAYER EQUATIONS IN PLANE OF 

SYMMETRY 

BOUNDARY-LAYER EQUATIONS IN PLANE OF SYMMETRY 

In reference 7, it is shown that the dimensionless laminar- 
boundary-layer equations for supersonic flow over a circular 
cone are 



(b) Variation of total head with angle of attack. 

Figure 4. — Variation with angle of attack of total head measured in 
boundary layer of cone. M , 3.095; 0, 7.5°. 


[f+ie^s+fggj] /xx-^gx/x,+ f (ifx) 2 +2/xxx=0 (la) 

r , 1 v'(<p) , 2 "I 2 

[y^~30 V ^30 30 gxgXv 

h xfx ~id ^ +2sxxx=0 (lb) 

r+(/x) 2 +(^) 2 =ri+«i 2 +w! 2 (lc) 

2>=^ f>T (Id) 

Equations (la) and (lb) are momentum equations, equation 
(lc) is an energy balance, and equation (Id) is the equation 
of state. A complete list of symbols is provided in appendix 
A. The functions /(X, <p) and g(\,<p) are related to the two- 
component vector potential discussed in reference 6 and are 
defined according to the relations 

u=h) 

> ( 2 ) 
W=g\) 
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in a manner such as to satisfy the continuity equation 
identically. 

The coordinate X has been formed as follows: 




(3) 


The coordinate <p is the angle between the vertical plane of 
symmetry of the flow and any meridional plane of the body 
(fig. 1). Equations (1) and (3) imply that parabolic simi- 
larity of the Blasius type exists in meridional planes. As 
pointed out in references 5 and 6, this conclusion of parabolic 
similarity applies for the boundary layer on any smooth 
conically symmetric body in supersonic flow (for example, a 
cone of elliptic cross section). 

In reference 7, all quantities are made dimensionless by 
referring them to the properties of the non viscous flow at 
the outer edge of the boundary layer when the cone is at 
zero angle of attack. In the present report it will be con- 
venient to use a different reference condition (subscript r ) 
which will be defined subsequently. The following quan- 
tities on the left are to be identified with the dimensionless 
groups on the right: 


u ) W'^'ulu T ,wlu T 
p^p/pt 
T ~ 2c p Tju r 2 


(4) 


velocity components should take on the corresponding 
nonviscous values 

/x(°°, *>)==%(#>) (6a) 

g\( m ,<p)=Wi(<p) (6b) 

At the cone surface, the u and w velocities should vanish 

/a ( o ,<p) ~g\(0j<p)=0 (6c) 

and the normal velocity v should vanish. It is shown in 
reference 6 that this last requirement is met if 

/(0,^)=g(0,p)=0 (6d) 

Equations (1), involving two independent variables, would 
be quite difficult to solve in general. However, a certain 
amount of information can be obtained by restricting con- 
sideration to the plane of symmetry, thus yielding a tractable 
set of ordinary equations involving X as the only independent 
variable. 

In the plane of symmetry (<p= 0,7r), io—g\— 0. Because 
u=j\ is even about the plane of symmetry and may be 
expected to be regular there, fx?— 0. The pressure and 
the density are also even, and therefore p' (<p) vanishes at 
<£>“0,7r. Thus, in the plane of symmetry, equation (la) 
reduces to the following equation: 

(V + 3^ /xx+2/xxx— 0 (7a) 


P^plp r Ur 2 

X^^PrUrZ/Cflr, p r U T y/Cfl r ^ 

where the constant C arises from the assumption of the 
temperature-viscosity relation of Chapman and Rubesin 
(ref. 9): 

Ih 0 ? ( 5a ) 

Mr -L r 


with C being defined as follows, in order to match equation 
(5a) to the, Sutherland formula at the cone surface (denoted 
by subscript w): 


r _fT w \i 1 +S/T r 
\Tj T w /T r + S/T T 


(5b) 


The quantity S may be taken as equal to (216° R)2 c P /u r 2 . 

The following additional physical assumptions are em- 
bodied in equations (1): 

(a) A thin boundary layer across which the static pressure 
is constant 

(b) Prandtl number of 1 and constant ratio of specific 
heats y 

(c) No heat transfer through the surface 

From equation (lc), since the case of Prandtl number of 1 
and no heat transfer is considered, T w in equation (5b) may 
be taken equal to the dimensionless stream stagnation 
temperature. 

The boundary conditions on the functions/(X,<p) and g(\,<p) 
are: At the outer edge of the bounda^ layer, the u and w 


Every term in equation (lb) vanishes at the plane of sym- 
metry; and, therefore, in order to obtain a meaningful 
equation, it is necessary first to differentiate equation (lb) 
with respect to <p and then drop terms which vanish at 
<p=0,7r. This procedure yields the following result: 

(/”^30 £^£xx*“^(gx*) 2 — £x*>/x— 20 ~ ~^+2£\xx*> = 0 (7b) 
Equation (lc) becomes 

T+(A) 2 =7 T 1 +U l 2 (7c) 

Equations (7) may be considered a set of ordinary differ- 
ential equations for the functions /(X,0) and g*(X, 0), or 
/(X,7r) and ^^(X,7r), depending on whether the solution is 
required at the bottom or the top of the cone. According 
to equations (2), the result /(X,0 or w) may be differentiated 
with respect to X to give the profile of meridional velocity u 
in the plane of symmetry. The form assumed by the cir- 
cumferential velocity profile w as <p -^> 0 or tt is given by 
£*>x(X, 0 or 7 r) in the sense that, at a small angular distance d(p 
away from the plane of symmetry, w ~g v \ d(p. 

The boundary conditions (eqs. (6)) become, in the plane 
of symmetry, 

/x(°o,0 or 7r)“Ui(0 or tt) \ 

^(“jO or 7 r)=«; 1 ^(0 or tt) > ( 8) 

/x(0,0 or ?r)=^x(0,0 or tt)=/( 0,0 or 7r)=^(0,0 or tt)=o) 

In view of the first of equations (8) and of equation (7c), it is 
convenient to specify the reference condition (subscript r) 
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to be that existing at the outer edge of the boundary layer, 
for the particular angle of attack under consideration, 
evaluated at either <p— 0 or <p— x, depending on whether the 
analysis pertains to the bottom or top of the cone. 

Because the pressure is assumed constant across the 
boundary layer, equations (Id) and (7c) and the assumption 
of constant pressure across the boundary layer (p—Pi) 
provide that 

i=^=l+^[l-(/x) 2 ] (9) 


From conventions (4), Ti= 



the following definitions are made: 


For convenience, 


£„(\,0 or tt)=y k\p(\) 


k=—w lv (0 or x) 


(10a) 

(10b) 


Equation (9) and definitions (10) are introduced into equa- 
tions (7a) and (7b), and a value of p"(0 or x) is assigned con- 
sistent with the nonviscous equations at the outer edge of 
the boundary layer (p"{0 or x) may conveniently be ob- 
tained from equations (lb) and (6) by setting g ipXX =g vX>iX =o 
when X= °o ). 

The following pair of simultaneous ordinary differential 
equations then results: 


(/+**)/"+2/'"=0 (11a) 

(f+wr+w"- w r 1 yr+ 

(*+|){ 1 +^r [1 - ( / /)2] } = ° (lib) 

and boundary conditions (8) become 

/ / (od) = ^'(oo) = 1 (12a) 

/'(0)=tf'(0) = 0 (12b) 

/( 0) = *(0)=0 (12c) 


Two parameters appear: k , which depends essentially on 
angle of attack, and Ti, which is essentially dependent on 
Mach number. If the angle of attack (and hence k) is zero 
or nearly zero, equations (11) become precisely those con- 
sidered in reference 7 and may be solved quite readily, since 
equation (lib) becomes linear and the solution of equation 
(11a) is well known as the Blasius function. When k differs 
substantially from zero (moderate or large angle of attack), 
equations (11) are both nonlinear and the solutions are 
interdependent. For any particular case, when only the 
stream Mach number, cone vertex angle, and angle of attack 
are specified, the parameters k and T\ must be obtained by 
recourse to a theory of the outer nonviscous flow. 

OUTER NONVISCOUS FLOW 

In references 1 to 3, the results of a theory of nonviscous 
supersonic flow about circular cones at angle of attack are 


tabulated. The case of zero angle of attack (ref. 1) is 
solved exactly in the sense that no assumption of small 
vertex angle is made. The equations are then expanded in 
powers of angle of attack with the use of the zero-angle-of- 
attack solution as the first approximation. Terms linear in 
angle of attack are presented in reference 2, and terms pro- 
portional to the square of angle of attack are presented in 
reference 3. Reference 4 clarifies the application of the 
theory to the computation of flow conditions at the cone 
surface. 

There are two objections to the use of this theory in the 
present application: 

(1) Neglect of terms in the expansion beyond that in- 
volving the square of angle of attack may lead to an insuffi- 
ciently accurate representation of the flow at the large angles 
of attack which are of interest. Unfortunately, no non- 
viscous theory is available that treats the effect of angle of 
attack with greater precision. In reference 4, a comparison 
of the theoretical and experimental pressure distributions is 
presented for a cone of semivertex angle of 10°, at a Mach 
number of 2, and of an angle of attack of 12.2°. The agree- 
ment shown is very good, especially since the angle of attack 
is sufficiently large that the pressure distribution is of the 
type shown in figure 2(b). 

(2) In reference 10, Ferri points out that the method of 
expansion used in references 2 and 3 is improper near the 
cone surface and leads to an erroneous form of the entropy 
distribution around the cone. Therefore, the theory cannot 
be applied if the vertex angle, the angle of attack, or the 
Mach number is so large that the flow may not be considered 
essentially isentropic. In reference 7, an argument is pre- 
sented to the effect that in the limit of infinitesimal angle of 
attack the presence of a boundary layer ensures that the 
error in entropy distribution is of no consequence even for 
large cone vertex angles. That argument in no way applies 
to the present analysis because the angles of attack con- 
sidered are not infinitesimal. For the purposes of this 
report, the use of references 2 and 3 in their present form is 
justified only in cases for which isentropic flow may be 
assumed. 

According to reference 4, the velocity components at the 
cone surface are, using the notation of references 1 to 3 for 
quantities tabulated therein, 


X 

— =^ = l + a = cos (p~\- 
u u 


a l? + ^ 5 cot e+ i + (H I cot e+ £) cos 2 *]+ • • • 


(13) 


w i(*>)_ 


u 


= a = sin <p+a 2 — csc0 — *== cot 0 sin 2<p-\- . . . (14) 
u ^ u u J 


The pressure and density are 

1 +* l cos *+« 2 [©+? M*+\l cot e+ 

V V Lv 2 2 p 

(f +5 w ~\ I cot e ) cos 2 ^]+ • • • (l5) 
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1 + a 4 cos v>+ c? [" M 2 +\ 4 cot 0+ 

P p L P * z P 

{fj+\ ^ 2 ~\ | cot e ) cos 2 ^J+ • • • ( 16 ) 

The barred quantities are those pertaining to the case of 
zero angle of attack. From equations (10b), (13), and (14), 


or, approximately, 


j 2 r , «2.o 2 /w 2 

k =zl ± ez +2 « (.W 


i yr^< 


£ 1 3g\~ 

m 26 w)_ 


+ 


(17) 


* 




2 ±a J +2a2 ( 

'w 2 

1 Vl-0 2 2\ 


d 6 u) . 

30 , , x. 

« 2 I 

,«o, u 2 \ ' 

\ u^~u) 

1 M 2 u 2 / 


2 — — 2 [”l i a — “h °? ( 

V 2 p pi_ u 2 L u \ 

TM 2 1-« 2 V 


'Pi p 1 —u 2 


The plus or minus sign refers to y = 0 or ir, respectively. 
From conventions (4), equations (13), (15), and (16), and 

the result of reference 1 that l/T=v?/(l— M 2 , 


. . 'M'O i «,MT, I I s/po . P2 l 2 « 2 \"1 

1 + v+v)l 1 1± “;+“ (? + 7 + ^i i=svJ 


+ 


1 ± 


»/ , 2 /P 0 , P2 , 27 S 2 \ 

a 3 +a yj+J+^iT^) 


or, approximately, 

1 W {, , /„*,{ , ,r o /l I «0 t ««\ , P0^P2 Po P2 o tt* 

r l 1 ±a l 2 5 + nj + ^L 2 v 1+ « + i j+T+r? f 2 i-i? +2 uVp pp_Ij + 


(18) 


In figure 5 are shown k and l/7\ as functions of a for a cone 
of semivertex angle of 7.5° and a stream Mach number of 
3.1. From the tabulations of reference 2, it may be inferred 
that under these conditions the isentropic assumption leads 
to errors of less than 1 percent in quantities proportional 
to the angle of attack. 

SOLUTION OF EQUATIONS AT *>=0 

Equations (11) have been solved, subject to boundary 
conditions (12), at <p=0, and various angles of attack for a 
cone of semivertex angle of 7.5° and a stream Mach number 
of 3.1 for which the values of k and 1/Ti are given in figure 5. 
The computations were carried out by Dr. Lynn Albers of the 
Lewis laboratory and are described in appendix B. The 
resulting boundary-layer profiles of meridional velocity u and 
gradient of circumferential velocity dw/d<p are shown in figure 
6. The curves for a=0 are obtained from reference 7. The 
profiles show clearly that, as the angle of attack is increased, 
the boundary layer becomes thinner on the bottom of the 
cone, and the shear stress at the wall increases. 

Skin friction. — The meridional and circumferential com- 
ponents of the viscous shear stress at the bottom of the cone 
surface may be written in coefficient form as follows: 




ipl-ttl 




where the quantities on the right are in dimensional form. 
Application of equations (Id), (2), (3), (4), (5), and (10a) 
yields 



M, 3.1; 0, 7.5°. 





(19b) 


y,.o=2/"(o) 


(19a) 


(19c) 
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Dimensionless coordinate, X 

Figure 6. — Velocity profiles at <p= 0. M } 3.1; 9, 7.5°. 


where 

£ — P lUlX 

Mi 

Variation of these skin-friction coefficients with angle of 
attack is shown in figure 7 for a particular case. 

Displacement thickness. — In reference 11 it is shown that 
the displacement thickness A for a cone at angle of attack is 
the solution of the equation 


iUi (A — 5^) + ^ [p 1 W 1 (A — 5^)] = 0 




At <p = 0, Wi = 0 and equation (20) may be solved directly 


dz~\-k8 <p 

1+k 


i 


( 22 ) 
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0 4 8 12 


Angle of attack, oi, deg 

Figure 7. — Meridional and circumferential skin-friction coefficients at 
<p—0. M, 3.1; 0, 7.5°. 



0 4 8 12 

Angle of attack, a, deg 

Figure 8. — Displacement thickness at y?=0. M , 3.1; 9, 7.5°. 


where k is defined in equation (10b). From equations (5) 
and (9), with 5 X defined in terms of a Reynolds number, 

r> PlUjlSa; 


[W'+^t (l -/ ,2) ] (23a) 
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Using equations- : {10) and applying FHospitaFs rule to 
evaluate at <p~ O’ |ihe limit of the ratio w/wi appearing in 
equation (21) yield | 

( 23 b) 

Therefore, 

yl(k X’lm <*+«]} 

(23c) 

Figure 8 shows th£ variation of displacement thickness with 
angle of attack for a particular case and again illustrates the 
expected progressive shift of the boundary layer from bottom 
to top as the angle of attack is decreased. 

LIMITATIONS OF METHOD AT <p=ir 

Except for quite small angles of attack, equations (11) 
cannot be solved at the top of the cone {ip— t). Over part 
of the range of angle of attack, the solutions are indetermi- 
nate; and, beyond a certain angle of attack, the solutions 
do not exist at all. These properties of equations (11) will 
be demonstrated and discussed in the following paragraphs. 

Asymptotic forms of equations. — The difficulties just 
mentioned may be$t be inferred from the asymptotic forms 
of equations (11) at large X. From equation (12a) it is clear 
that, for large X, jf and $ may be written as follows: 

f=\+F{\) 

^=X+T(X) 

where F' ( os ) = \f/'{ co ) = 0. Substitution into equations (11) 
yields the asymptotic forms for large X 

(1+A:)XF" + 2F"' = 0 (24a) 

(1 + 2 (fc+|) *'=| F' + A (* + §)] 

(24b) 


Indeterminate solutions. — Consideration will now be given 
to the problem of obtaining the complementary solution of 
equation (24b). Defining a new dependent variable 


yields the equation 

^"-[f fc+^+(^~ ) 2 X 2 ]^=° (25) 

This is essentially Weber’s equation (ref. 12, paragraph 
16.5), and the asymptotic solutions are 


e 



_3 fc+ 5 / 9 
i+fc 


e 



k± 1/3 
1-Hfc 


Thus, the asymptotic complementary solutions of equation 
(24b) are 


fe+5/9 


-¥”(/ t*x) ,+ * 

(26a) 

,'*+1/3 


■'■(V't* *) I+ ‘ 

(26b) 


Because it is required that 4>'(o°)— 0, solution (26b) is 
rejected when &>— ■ 1/3. If &< — 1/3, both solutions may be 
retained, and an additional undetermined constant appears. 
When 0 (at the bottom of the cone), the complete solu- 
tion of equations (11) exists and is unique, and solution 
(26b) is to be rejected in forming the asymptotic solution. 
Therefore, if solution (26b) and the associated constant must 
be retained when &< — 1/3, it is clear that the complete 
solution for — 1/3 cannot be unique. This indeterminacy 
has been verified by numerical integration of equations (11), 
as described in appendix B. 

This indeterminacy arises because essential information 
has been lost by specializing the equations to apply only in 
the plane of symmetry. When the equations are so written, 
it is implied that the boundary-layer development in the 
plane of symmetry is affected only by conditions in that 
plane. The lateral region of influence of points on the body 
in the plane of symmetry grows parabolically (the shaded 
regions in sketches (a) and (b)), according to the law of 
molecular diffusion, when, as in the present instance, there 
is no pressure gradient in the stream direction. Fluid enters 
the boundary layer from the outer stream. If the fluid then 
moves laterally (because of the angle of attack) out of the 
region of influence of the plane of symmetry, as shown in 
sketch (a), the flow is uniquely determined by outer stream 
conditions in the plane of symmetry. Clearly, this is the 
case when <p—tt and a is small, and when <p = 0 and a has any 
positive value. When except for small angles of 

attack, the lateral motion of the fluid is inward relative to 
the region of influence, as shown in sketch (b). This fluid 
then brings into the region of influence of the plane of sym- 
metry information concerning boundary-layer development 
as it proceeds around the cone from the bottom. Conse- 
quently, outer stream conditions at <p—tt may not uniquely 
determine the boundary-layer characteristics at <p=7r, and 
indeterminate solutions of equations (11) may be anticipated. 



f 




f 

(c) 
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The condition k= — ){ specifies an angle of attack such that, 
at <p= t, the outer flow streamlines are just tangent to the 
parabolic region of influence. This statement may be 
proved as follows: From definition (10b), if k— — % , then 
WiJui=—d/2. Near w (see sketch (c)), Wi=Wi p d<p= 

Wi —• Thus, when k=—%, Wi/Ui~t/2. From this last 

equation and the geometrical properties of parabolas, it may 
be inferred that the streamlines are parabolas with focus at 
the cone apex. Therefore, the situation shown in sketch 
(b) applies if k<^ — ){ and equations (11) have no unique 
solution. 

It is noteworthy that the mathematical nature of the 
indeterminacy so far described is quite similar to that which 
characterizes solutions of the Falkner-Skan equation for the 
plane boundary layer when the outer flow velocity is pro- 
portional to x~ m ( x being the streamwise coordinate and m a 
positive number). Hartree (ref. 13) has treated this prob- 
lem in detail and finds an “extra” asymptotic solution similar 
to equation (26b), decreasing with distance away from the 
surface according to a power law. This extra solution is 
rejected by Hartree on the ground that an unrealistically 
thick boundary layer would be predicted if the solution 
were retained. 

A thorough study of the asymptotic nature of equations 
(1) might lead to a similar conclusion in the present problem. 
However, in the present problem, the mathematical difficulty 
goes deeper; when k^> — %, numerical integration failed to 
give unique results, even though solution (26b) is certainly 
to be rejected. Thus, the condition k> — % is not sufficient, 
but rather only necessary for uniqueness. 

Inspection of the profiles shown in figure 6 for a=k—0 indi- 
cates that, for 0>&>— %, the streamlines within the bound- 
ary layer may be expected to incline more sharply toward 
the plane (p= tt than do the streamlines near the edge of the 
boundary layer and thus may bring information from beneath 
the cone even though the outer ones do not. Therefore, the 
necessary condition for uniqueness would be (see eqs. (2) 
and (10a)): 



Figure 6 indicates that perhaps the maximum value of If 
is to be found at X=0, in which case criterion (27a) would 
become 


Nonexistence. — The solution of equation (24a) is 


l+fc 

F”=e 4 




The requirement that F"(a °) = 0 is met only if k^> — 1. If 
k< — 1, no solution of equation (24a) exists which satisfies 


the boundary conditions, and, therefore, the Prandtl boundary- 
layer equations fail to describe the flow. This was first 
pointed out by Hayes in reference 5. The Prandtl equations 
differ from the exact equations essentially in that a thin 
boundary layer is assumed. Thus, if k<~ 1, the boundary 
layer cannot be regarded as thin. It may be noted that 
equation (22) implies that, as Ar — > — 1, the displacement 
thickness approaches infinity. 

Any boundary layer grows by the entrainment of fluid at 
its o\iter edge. That is, fluid particles acquire vorticity by 
entry into the boundary layer. The reverse process cannot 
occur — fluid particles cannot leave the boundary layer, thus 
losing their vorticity. In the case under consideration, it 
will be shown that when — 1, the streamlines at the outer 
edge of the boundary layer would proceed outward relative 
to the boundary layer, if the boundary layer were to retain 
parabolic similarity. Because such a situation is physically 
impossible, the Prandtl equations fail to yield a solution. 

From reference 11, the normal velocity at the outer edge 
of the boundary layer at <p = tt is 

»,=<«, £j+(*-a> (!■•?),_„ (28) 

where A is a somewhat arbitrary definition of the outer edge 

of the boundary layer and (^-^) is obtained from analysis 

of the outer nonviscous flow. The equation of continuity 
for the outer flow, evaluated at the surface of the cone in the 
plane of symmetry, may be written 


_l 

L J_ 

bw ' i ( 

X 

r dx 

1 \ 


(tL=° 


(29) 


With equations (28) and (29) combined, the flow inclination 
at the outer edge of the boundary layer is 


Vi_dA h—A / , 1_ dwA 

Ui~~ dx x V ~^dui d<p J 


(30) 


The Prandtl (thin) boundary layer may be expected to exist 
only if 

bh.Vi 
bar Ui 


or, with equation (30) introduced, if 

If parabolic similarity is assumed ( h and A each proportional 
to fx), inequality (31) becomes 

1 dw x ^ 3 

6u j d<p ^ 2 

or, from equation (10b), 

k^> — 1 


ii i i 


Mill 


llll 
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SEPARATION 

The critical condition A= — 1 might be expected to be of 
physical as well as analytical significance because some sort 
of catastrophic thickening of the boundary layer is implied. 
In particular, this critical condition may reasonably be 
supposed to be connected with the phenomenon of laminar 
separation. In order to explore this possibility, it is first 
necessary to describe qualitatively what is meant by separa- 
tion in three-dimensional boundary-layer flow. Difficulty 
has been encountered in establishing a satisfactory qualitative 
criterion for three-dimensional separation (see, for example, 
ref. 5) . Therefore, in the subsequent paragraphs, the 
general problem of three-dimensional separation will be 
discussed, and then the particular case of the cone at angle 
of attack will be considered. 

GENERAL CONSIDERATIONS 

Plane flow. — In plane flow, separation is customarily 
identified by the appearance of reverse flow (sketch (d)). 
In order to generalize this concept to three-dimensional 
flow, it is necessary to consider the separated region as a 
whole. In plane flow, the separation point of sketch (d) 
might be regarded as the forward boundary of a vortex 
sheet embedded, or encapsulated, within a region bounded 
by the body and a stream surface meeting the body (sketch 
(e)). Sketch (e) shows separation followed by reattach- 
ment. Of course, the sort of separation of greatest engineer- 
ing importance occurs when such an embedded vortex sheet 
rolls up to form a large concentrated vortex, or is shed as a 
vortex street, with the consequence that the outer flow is 
greatly disturbed and a large pressure effect (form drag) 
occurs. 




If a boundary-layer solution of the type shown in sketch 
(e) were obtained, the vortex sheet would be completely 
embedded in a thin boundary layer, and would presumably 
tend to remain flat against the body (leaving out of account 
possible effects of laminar instability). However, if the 
theoretically predicted vortex sheet extends downstream 
into a region where the thin boundary-layer equations do 
not apply (that is, where the solution “blows up” predicting 
an infinitely thick boundary layer), then in that region the 
vortex sheet would not be constrained to lie flat and the 
rolling-up process would occur. 

The foregoing discussion seems to provide the proper 
description of what occurs behind a bluff body: A complete 
solution of the thin laminar boundary-layer equations for 
the flow over a cylinder would probably yield a streamline 



pattern of the form shown in sketch (f). The boundary 
layer would be predicted to gain mass flow by entrainment 
as it proceeds around the body, until it reaches the vicinity 
of the rear stagnation point. There, the mass flow contained 
in the boundary layer must finally leave the body and proceed 
downstream. The boundary layer therefore cannot remain 
theoretically thin, but rather must approach infinite thick- 
ness in violation of the Prandtl assumptions. In this 
region, then, the aft boundaries of the pair of vortex sheets 
are free to roll up into concentrated vortices, thus distorting 
the outer flow in such a way that the rolling-up process 
engulfs most of the region which would otherwise be occupied 
by a flat vortex sheet. The leading edge of the sheet, 
however, is still constrained to lie flat against the body. 

Three-dimensional flow. — The foregoing description of 
plane separation may be generalized to three-dimensional 
flow as follows: A separated region on a three-dimensional 
body consists of a vortex sheet embedded between the body 
surface and a stream surface attached to the body in a closed 
curve, as shown in sketch (g), which is a view of the body 
from above. The arrows indicate possible directions of 
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resultant shear stress at the surface and outside the separated 
region. The situation shown in the sketch would correspond 
to separation and reattachment in plane flow. If somewhere 
within or at the boundary of the separated region the theo- 
retical boundary-layer solution would a priori be expected 
to blow up, then the vortex system within the separated 
region is free to roll up into a more or less vigorous system of 
vortices. 

Thus, thin-boundary-layer theory can be used to obtain 
the following information concerning laminar separation: 

(1) The solution may establish the existence of a vortex 
sheet which is embedded in a flat bubble on the surface and 
which could adhere to the surface and remain part of a thin 
boundary layer, provided that the Prandtl equations are 
valid everywhere in the separated region. 

(2) The solution might predict the boundary layer to go 
to infinite thickness somewhere in the separated region. If 
this happens, then the separated region is free to roll up, thus 
providing a vigorous wake (which, of course, is not amenable 
to boundary-layer theory). 

For flow about a plane body, the boundary-layer solution 
is not needed for predicting the breakdown of the Prandtl 
assumptions. Physical considerations suffice to establish 
where (at stagnation point of outer flow) and when (always) 
the breakdown occurs. In three-dimensional flow this is not 
always so clear, and at least certain features of the solution 
are required to be known. In order that the boundary-layer 
equations be applicable, the solution must be such that the 
boundary layer entrains fluid (that is, flow streamlines enter, 
but do not leave, the boundary layer at its outer edge). In 
a Cartesian system, where h(x,z) is the outer edge of the 
boundary layer, this requirement may be written as 


Uih x -^Wihz^>{v\) y =h (32) 

Equations (32) and (28) may be combined with the equation 
of continuity in the form 

Pi = — [(pi^l)*+(pi^l)z] 

to yield 

A p-a^i+A p-A) Pl w 1 ]>o 


or, in vector notation, 

div [(A—A)pigi]>0 (33) 

where A is the displacement thickness, and is the velocity 
vector in the outer flow evaluated at the body surface. 


In many cases, circumstances may be found for which 
inequality (33) cannot be satisfied. For example, for plane 
incompressible flow about a cylinder and, as is customary, 
with u and x defined parallel to the surface, inequality (33) 
may be written 


dx 


(h~ A) 


> 


du\ 

dx 

Ui 


As the rear stagnation point of the outer flow is approached, 
Ui tends to zero, while —dujdx remains finite and positive. 
Therefore, since h — A by definition must be greater than 
zero, d(h—A)/dx must approach infinity in clear violation 
of the boundary-layer assumptions. 

SEPARATION ON CONE AT ANGLE OF ATTACK 

In the previous discussion, it was concluded that separation 
involving a strong vortex pattern occurs if a tentative 
boundary-layer solution predicts an embedded vortex sheet 
coupled with a local breakdown of the assumption of a thin 
boundary layer. In the case of the cone, inequality (33) 
may be used to predict the circumstances under which the 
boundary layer may not be regarded as thin: When the fact 
of parabolic similarity is introduced (h — A proportional 
to --yfx), inequality (33) becomes 




h — A 


Wi 


(34) 


This inequality indicates infinite boundary-layer thick- 
ness only when W\ — 0 and, even then, only if W\<p is negative 
(which is true at the top of the cone, <p= ir) and larger in 
magnitude than \Qu\. This is true only for angles of attack 
larger than that for which ft— — 1 (by eq. (10b)). When 
the angle of attack is smaller than this critical value, the 
right member of inequality (34) is always negative, and 
d(A— A) [dtp may be considered to vanish by symmetry 
at <p=ir without violating inequality (34). 

The foregoing result may be explained on physical grounds 
as follows : As the boundary layer proceeds around the cone, 
it entrains fluid which it then conveys toward the top 
(symmetrically, from both sides of the cone). In the plane 
cylinder case, the fluid similarly conveyed must finally 
erupt from the boundary layer when the stagnation point 
is reached. However, on the cone, the boundary layer grows 
parabolically along generators; and, hence, if the crossflow 
is not too strong (small angle of attack), the fluid brought 
to the top may simply become part of the growing boundary 
layer. For larger angles of attack, the boundary layer 
cannot grow at a rate sufficient to absorb the additional 
fluid, and eruption occurs with the consequent breakdown 
of the thin-boundary-layer assumptions. 

Accordingly, it is proposed that when the angle of attack 
is less than that for which k= — 1, a thin boundary layer 
may cover the cone (fig. 3(a)). For larger angles of attack, 
any vortex sheet present will roll up to form attached lobes 
(fig. 3(b)); for still larger angles of attack, a vortex street 
is produced. 

Thus, when — 1 (angle of attack greater than that for 
which k= 1), strong viscous cross forces (viscous lift) on the 
cone may be expected. These forces are discussed by Allen 
and Perkins in reference 14. Of course, a weakness of the 
present analysis is that no indication is given as to the 
strength of the rolled-up vortex system because, when k— — 1 , 
the presence of an embedded vortex sheet over the top part 
of the cone has not been established. It seems likely that 
such a vortex sheet does exist because a rather strong adverse 


h — A 
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Figure 9. — Minimum angle of attack for which separation appears at 


<P~7T. 

pressure gradient (fig. 2(b)) always exists when k— — 1. In 
fact, it may be shown (most conveniently by evaluating eq. 
(7 a) at \= oo ) that an adverse pressure gradient exists 
when &< — 2/3. 

Equations (13), (14), and (10b) give the critical angle of 
attack corresponding to &= — 1. This angle is presented as 
a function of stream Mach number and vertex angle. Figure 
9 shows the results of such a calculation. The critical angle 
of attack is given as a ratio of angle of attack to semivertex 
angle for convenience. The results suggest that, in general, 
separation involving lobes occurs later (in terms of relative 
angle of attack a/0) for the smaller vertex angles, particu- 
larly at higher Mach numbers. Figure 9 indicates the pos- 
sibility of rather profound qualitative differences in the 
flows at high Mach number about cones of different vertex 
angles. 

The foregoing interpretation of the critical condition 
k= — 1 is supported by the experimental result shown in 
figure 4. From figure 5, Jc= — 1 when a— 6.2°, under the 
conditions of the test. Figure 4(b) shows the measured 
total head rising as the angle of attack is increased beyond 
6.2°. Possibly this effect is caused by the induced field of 
the symmetrical pair of vortex lobes sweeping away the 
thick boundary layer between, thus reestablishing a thin 
boundary layer at the top of the cone. 

It may be of interest to note that if similarity also holds 
for the turbulent boundary layer on a cone, and the similarity 
law is nearly linear (rather than parabolic as in the laminar 
case), separation would first appear at a higher angle of 
attack than in laminar flow. In fact, equations (10b) and 
(31) or (33) would yield the criterion Jc= — 4/3. 

For the boundary layer produced by supersonic flow over 
any smooth conically symmetric body in supersonic flow 
(such as a cone of elliptic cross section), inequality (33) and 
the condition of parabolic boundary-layer similarity may be 
used to find a criterion equivalent to k= ~ 1 for the maximum 
angle of attack consistent with a thin boundary layer. 


CONCLUSIONS 

The laminar boundary-layer flow about a circular cone at 
large angles of attack to a supersonic stream has been ana- 
lyzed in the plane of symmetry with the following results: 

1. At the bottom of the cone, profiles of meridional ve- 
locity and of the gradient of circumferential velocity were 
determined and showed the expected tendency of the bound- 
ary layer to become thinner on the underside of the cone as 
the angle of attack is increased. 

2. At the top of the cone, except for very small angles of 
attack, the analysis (which is restricted to the plane of sym- 
metry) failed for the following reasons: 

(a) For angles of attack greater than some rather small 
value, the boundary layer brings information from beneath 
the cone into the vicinity of the plane of symmetry at the top. 
Therefore, the analysis, which deals only with the plane of 
symmetry, yielded indeterminate solutions. 

(b) For angles of attack greater than some angle (roughly 
of the order of the cone semivertex angle), no boundary- 
layer solution is possible. The characteristics of the outer 
flow and the known parabolic similarity of the boundary 
layer would together imply that, beyond this critical angle, 
there would be a component of flow leaving the boundary 
layer. This is physically impossible, since a boundary layer 
always entrains fluid. Thus, beyond the critical angle of 
attack, no solution can exist for equations which presume a 
thin boundary layer. 

For three-dimensional flow it is proposed that a separated 
region be regarded as a vortex sheet embedded in the bound- 
ary layer, remaining flat against the body if the assumption 
of a thin boundary layer is valid throughout the region. If, 
however, the boundary-layer assumptions break down any- 
where in the separated region, it is inferred that the vortex 
sheet may roll up to form strong vortices which may either 
remain attached or be shed as a vortex street. 

On the cone, therefore, the critical angle of attack beyond 
which no boundary -layer solution is possible at the top of 
the cone represents the maximum angle of attack for which 
the boundary layer is everywhere thin or, alternatively, the 
minimum angle of attack for which major disruption of the 
flow may be expected because of the formation of strong 
vortex lobes. Beyond this angle of attack, strong viscous 
cross forces may be anticipated. 

A similar criterion could easily be obtained for the bound- 
ary layer on any smooth conically symmetric body in super- 
sonic flow. 

The assumption of a suitable similarity law suffices to es- 
tablish a similar criterion if the boundary layer on a conical 
body is turbulent. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 15 , 1952 
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APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 


c 

constant appearing in temperature-viscosity 
relation (eq. (5a)) 

Of. 

component of skin-friction coefficient in 
x-direction 

Of. 

component of skin-friction coefficient in 
^-direction 

Cp 

specific heat at constant pressure 

F(k) 

function appearing in asymptotic representa- 
tion of / (eqs. (24)) 

fM 

function related to meridional velocity u by 
eq. (2) 

g(\<p) 

function related to circumferential velocity 
w by eq. (2) 

h 

height of outer edge of boundary layer 

k 

related to circumferential gradient of cir- 
cumferential velocity in plane of symmetry 
(eq. (10b)) 

M 

Mach number 

V 

static pressure 

gi 

velocity vector at outer edge of boundary 

layer 

R x 

Reynolds number, piUix/pi 

R&j R& x y Rip 

Reynolds numbers, p^A/jUi, Pi^xhlPu 
PiUidp/pi, respectively 

T 

absolute static temperature 

u 

meridional component of velocity 

V 

component of velocity normal to surface 

w 

circumferential velocity component 

X 

coordinate along generators of cone 

y 

coordinate normal to surface 

a 

angle of attack (positive as shown in fig. 1) 


y 

A 

5 * 

5* 

0 

o 

x 

M 

P 

<P 

*(X) 

*(X) 

Subscripts: 

max 

r 


1 


Superscripts: 


ratio of specific heats 
displacement surface height 
mass-flow defect associated with meridional 
velocity profile (eq. (21)) 
mass-flow defect associated with circum- 
ferential profile (eq. (21)) 
semivertex angle of cone 
sine of semivertex angle of cone 
dimensionless variable (eq. (3)) 
coefficient of viscosity 
density 

angular coordinate around cone 
function appearing in asymptotic representa- 
tion of ^ (eqs. (24)) 

function related to circumferential velocity 
w in plane of symmetry by equation (10a) 

maximum 

reference condition, nonviscous flow at 
surface, at <p—0 or <p— 7r, whichever is 
appropriate 

evaluation at outer edge of boundary layer 
(alternatively, nonviscous flow at surface) 
Subscript notation for partial differentiation 
has been used 

Primes denote ordinary differentiation with 
respect to X or <p 

Bar over quantity indicates evaluation of 
nonviscous flow at surface when cone is 
at zero angle of attack 


APPENDIX B 

NUMERICAL SOLUTION OF DIFFERENTIAL EQUATIONS 

By Lynn Albers 


The two simultaneous nonlinear ordinary differential 
equations (11a) and (lib) together with boundary conditions 
(12) constitute a two-point boundary-value problem. The 
method of numerical solution used applies directly only to 
problems for which all boundary conditions are specified at 
a single initial point (the origin, in the present case). Each 
numerical integration was therefore performed starting with 
boundary conditions (12a) and (12b) and a tentative speci- 
fication of/"(0) and ^"(0). In each case, such integration 
was carried out for a sufficient variety of conditions /" (0) 
and ^"(0) so that the correct set of initial conditions yielding 
the proper behavior at X= oo (boundary condition (12a)) 
could be inferred to the desired degree of accuracy. 

Integration was performed according to the following basic 
scheme: With the value of /"'(X) and ^'"(X) given at five 
closely spaced values of X, fourth-degree polynomials may be 


passed through the two sets of values off" and Then, 

if f } f } /", \f/ } \j/', and \p" are known at the fifth point, the 
polynomial representations of f" and may be integrated 
to yield /, /', /", and \p" at the next (sixth) point. 

These quantities may then be substituted into differential 
equations (11) to yield/"' and \p'" at the sixth point. In 
this way, the solution may be extended one step at a time, in 
each step by use of the solution at the five previous points. 
In order to begin this procedure, the solution must first be 
found at five points starting at the origin and must be subject 
to boundary conditions (12b) and (12c) and the tentative 
selection of/"(0) and ^"(0). 

This preliminary calculation was done in the following 
manner: /'"( 0) was calculated directly from equations (11) 
and (12) and was used as an initial estimate of/'" at the next 
four points. Given /(0),/'(0), and/"(0), the values of/,/', 
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and/" were computed at the second point by integrating a 
fourth-degree polynomial passed through the five values of 
/'". In a similar manner, \f/ f \f/', and at the second point 
were found. Direct substitution into equations (11) then 
yields improved estimates for/"' and \p"' at the second point 
and thus an improved polynomial representation of these 
functions which may be used to obtain values of/,/',/", 

V > and at the second point, and so forth, until improved 
values have been obtained at the fifth point. This procedure 
was repeated in an iterative manner until convergence was 
obtained at each of the five initial points. 

All calculations were performed on the IBM Card Pro- 
grammed Electronic Calculator. Results are considered 
correct to four significant figures. 
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